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Abstract 



A spin system is considered with a Hamiltonian typical of molecular magnets, 
having dipole-dipole interactions and a single-site magnetic anisotropy. In addition, 
spin interactions through the common radiation field are included. A fully quantum- 
mechanical derivation of the collective radiation rate is presented. An effective 
narrowing of the dipole-dipole attenuation, due to high spin polarization is taken 
into account. The influence of the radiation rate on spin dynamics is carefully 
analysed. It is shown that this influence is completely negligible. No noticeable 
collective effects, such as superradiance, can appear in molecular magnets, being 
caused by electromagnetic spin radiation. Spin superradiance can arise in molecular 
magnets only when these are coupled to a resonant electric circuit, as has been 
suggested earlier by one of the authors in Laser Phys. 12, 1089 (2002). 
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1 Introduction 



Superradiance is a phenomenon well known for atomic systems [1,2], where it occurs be- 
cause of the self-organized correlation of atomic transitions through the common radiation 
field. There exists a similar phenomenon of spin superradiance arising in spin systems, 
as is reviewed in Refs. [3,4]. An accurate microscopic theory of spin superradiance has 
recently been developed [5-11] demonstrating that, despite many similarities between the 
atomic and spin superradiance, these phenomena possess features that radically distin- 
guish the spin from atomic superradiance [3-11]. The main such a fundamental difference 
is that the spin superradiance in real systems cannot be caused by their magnetodipole 
radiation field. The fact that the spin radiation rate is negligible as compared to all other 
relaxation rates was, first, noticed by Bloembergen [12]. This radiation rate cannot yield 
the development of spin superradiance, as has been emphasized recently [13] and analysed 
in detail in review [4]. In order to realize spin superradiance, it is necessary to couple 
a polarized magnet to a resonant electric circuit. The collectivization of spin motion 
happens due to the resonator feedback field, which, thus, replaces the common radiation 
field [3-11]. The analogies and differences in atomic and spin superradiance have been 
extensively analysed in reviews [3,4,14]. 

Even though the physics of spin superradiance has been well established, there have 
recently appeared discussions on the role of the radiation rate in spin relaxation. The 
interest to this problem has been resumed in connection to the studies of a novel class 
of materials, called molecular magnets, which can be highly polarized and may possess 
large molecular spins (see review [4]). The possibility of realizing spin superradiance by 
molecular magnets, coupled to a resonator, was advanced in Ref. [10], being based on 
the well developed theory [3-11]. At the same time, there have appeared speculations 
that spin superradiance in molecular magnets could be produced without coupling them 
to a resonant circuit, but in a way completely equivalent to atomic superradiance, so 
that the collectivization in spin motion would be due solely to their interactions through 
the common radiation field. Plausible collective effects in molecular magnets, caused by 
the radiation interaction, were discussed in Ref. [15], though the consideration there 
involved a very strong transverse magnetic field. At this instance, it is necessary to stress 
that superradiance as such is defined as a self-organized process, without being pushed by 
strong transverse fields. This definition was given in the original paper by Dicke [16] and is 
generally accepted for both atomic [1,2,17] as well as spin systems [3,4]. Having a generally 
accepted definition, there is no reason to change it. Therefore a process, involving strong 
transverse fields, has nothing to do with superradiance. In the best case, this could 
be collective spin induction, provided there would really be any collective effects. The 
possibility of the collective spin induction, caused by the molecular interactions through 
the electromagnetic radiation they are emitting, was considered in Ref. [15] in the frame of 
the phenomenological Landau-Lifshits-type equation for the classical magnetic moment, 
without taking account of the dipole spin interactions. However these dipole interactions 
do exist and are very strong in all really available molecular magnets. 

An effort of taking into account the dipole interactions was attempted in Ref. [18]. 
Unfortunately, the authors confused the effective spin temperature with the real tempera- 
ture of the lattice. These temperatures as is well known [19-21] have nothing in common, 
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when the system is not in equilibrium. The lattice temperature is intrinsically positive, 
while the spin temperature is positive or negative depending upon the average expectation 
value of spin polarization. Moreover, there even exist two effective spin temperatures for 
the same spin system, the Zeeman and dipole effective temperatures. These temperatures 
are directly related to the average spin, which in a nonequilibrium system is a function 
of time. As a consequence, the effective spin temperatures are also functions of time and 
are defined by the Provotorov evolution equations. For inverted nonequilibrium spins, the 
Zeeman temperature is negative. One should not confuse the negative time-dependent 
effective Zeeman spin temperature with the positive stationary real lattice temperature. 

In the present paper, we give an accurate analysis of the dynamics for a spin sys- 
tem described by a microscopic Hamiltonian typical of molecular magnets, including the 
dipole interactions together with the electromagnetic-field interaction. All consideration is 
based on a fully quantum-mechanical picture and on well grounded mathematical meth- 
ods of treating nonlinear spin-evolution equations. We do not couple the magnet to a 
resonant circuit, which allows us to concentrate our attention on the role of the radiation 
interactions and on the influence of the related radiation rate on spin motion. Our firm 
conclusion is that these radiation interactions play no role and the radiation rate is negli- 
gible, but the spin dynamics is governed by the dipole spin interactions. In this way, the 
spin interaction through the radiation field is absolutely unable to produce superradiance 
in molecular magnets. 



2 Electromagnetic spin interactions 

A realistic microscopic Hamiltonian for molecular magnets has the form 

N 1 JV 

£ = Effi + ? E£y + £/. (1) 

i=i z i+j 

Here the first term describes the single-spin energy, with 

Hi = -fioS* • B, - D{Stf , (2) 

where /xo = ^7s, 75 being the gyromagnetic ratio for a molecule of spin S; D is a single-site 
anisotropy parameter; the total magnetic field 

B, = B e z + H, , (3) 

acting on each spin, contains an external magnetic field Bq and the radiation field Hj = 
H(rj,t). The second term in Eq. (1) corresponds to dipole spin interactions, with 

^ = E<^f, (4) 

a/3 

where 

' if' 
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is the dipolar tensor and 



II 

The last term of Eq. (1) is the Hamiltonian of electromagnetic field 

^ = i-/(E 2 + H 2 ) dr, (5) 

with an electric field E = E(r,t) and magnetic radiation field H = H(r, t). The latter is 
expressed through the vector potential A as H = V x A. In what follows, the Coulomb 
calibration V • A = is used. 

To write the equations of motion for the spin operators Sj in a compact form, it is 
convenient to introduce the notation for the local fluctuating fields 



£o = \ E + + cyst) , 



h 

^ = £ E (2cy5/ - ^ OyfiJ" + 26yS+) , (6) 

in which Sj* 3 are the ladder spin operators and 

n ..= r)zz h - — - ( D xx - D vy - 2iD x A r ■ — - ( D xz - iD yz \ 

The magnetic moment operator can be represented as 

Mi = /2 Si = flSt + F s i + VoS* , (7) 
with the employed notation 

\i = — (e x - ie y ) , fj, Q = fi e z . 

The latter vectors enjoy the properties 

l/2| 2 = f, H 2 = 0, /2-/2 = 0. 
And two more important notations are the Zeeman frequency 

ujq = - i n B (8) 

and the effective field 

/=_ * + (9) 

Then the equations of motion for the spin operators read as 
dSi 



dt 



-i(wb- jrik-Hi + t )sr + fS? + j:D(SrS? + S?Sr) , 
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dSf 



- \ {f + s- + Sff) 



dt 2 

The equations for the electromagnetic field are 

1 <9E -* Aty 1 OA 

- — = VxH j, __ = _e 

c dt c c dt 

where the current density is 

N 

j = -c^M,xVi(r-r'). 

i=i 

From Eqs. (11) one has the equation 

(v 



(10) 



'ID 



(12) 



3 i^A = "^J, 



whose solution is 



c 2 <% 2 ; 



A(r,i) = - / j r',t 



r — r 



dr 7 



|r — r'l 



(13) 



(14) 



The average vacuum fluctuations are assumed to be zero. The vector potential (14), with 
the current density (12), can be represented as the sum 



in which 



Aj — Aj + A.f + Aj , 



1 d\ r, 



(15) 



c dt / rh 3 



AT = -E i + 

4 = -E£x/*,^(t- ^) , 



(16) 



where the property 



d 



dr 1 



S~ t- 



c dt b * V c 



is used. From Eqs. (15) and (16), one gets 

H, = V, x A, = H- + H+ + , 

with the expressions 



(17) 



HT = -E 



fl* - (g* ■ n^)n^- <9 2 /2* - 3(/2* ■ n^n^- / _d 



<r r. 



<9t 2 



c ^ 



c 



(18) 
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Keeping in mind a macroscopic in all directions sample, we may simplify Eqs. (18) 
by averaging them over the spherical angle f2( n «j) related to the unit vector n^-. In doing 
this, we employ the integrals 



If 2 

— j[fl-(jl- n)n] dQ(n) = - fl , 

— f[jl- 3(/2- n)n] dft(n) = , 

47T 7 

If 1 

— y (//• n)n dfi(n) = - p, . 



4tt 

Then from Eqs. (18), we find 

^ = -\FY.~s-(t-^) (19) 

and = 0. 

Since we are looking for the possibility of arising coherent collective effects, we have to 
accept the first basic condition for the existence of such effects. This necessary condition 
is the assumption that the transverse spin motion can be characterized by a well defined 
frequency u, with the related wavelength A = 2n/k and wave vector k — uj/c. Under this 
condition, the Born approximation is valid, 

Sj(t- ^=Sr(t)&(ct-r)e lkr , (20) 
where O(-) is a unit-step function, which allows us to rewrite Eq. (19) in the form 

H" = \ feV E e ^T^ 9(<* - r l3 )Sj{t) . (21) 

3 % J 

The magnetic field (17), with taking account of Eq. (21) and of the properties jl • = 
and /2o-H~ = 0, has to be substituted to the effective force (9) entering the spin equations 
of motion (10). The latter equations can be treated by the scale separation approach 
[5-10,14,22], whose mathematical foundation is based on the averaging techniques [23]. 
In this approach, one treats the local fluctuating fields (6) as random variables, which 
allows for the account of quantum spin fluctuations due to dipole interactions. Then the 
equations of motion (10) are to be averaged over the spin degrees of freedom. To this 
end, we define the averages for the transverse spin 

u(r,t) = ^ <S-(r,t)>, (22) 

coherence intensity 

w(r,t) = ± <S + (r,t)S-(r + 0,t)> , (23) 

and spin polarization 

s(r,t) = ± <S z (v,t)> . (24) 
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To proceed further, one needs to invoke the second fundamental condition for the 
occurrence of collective effects. This is the existence of the regions of strongly correlated 
spins, which are called spin packets. Really, without a substantial correlation, no collective 
effects can arise. The existence of a strongly correlated spin motion in some regions of 
the sample implies that the spin motion inside each of the regions is coherent, while the 
correlation between different regions of correlation is either absent or very weak. Let the 
characteristic length of a correlated region be L c , volume V c = L\, and the number of 
spins inside it N c . This characteristic length has to satisfy the inequality 

kL c < 1 . (25) 

If inequality (25) does not hold, then all spins in sum (21) oscillate independently of each 
other, with their correlation limited by only nearest neighbours, which follows from the 
fast spatial oscillations of the kernel in sum (21). But under condition (25), this kernel 
is practically constant, and all spins in the region V c can be synchronized. When the 
wavelength A ^> L is much larger than the system length, then L c = L, and the whole 
system radiates coherently. This, however, is not compulsory. The wavelength A can be 
much smaller than L. In that case, the system separates into several correlated regions, 
which radiate almost independently from each other. As a result, the total radiation pulse 
has an oscillatory behaviour. This is a well known picture in optical coherent radiation 
[1,2]. Lasers with a large aperture always radiate not by a unique beam but by a bunch of 
filaments [24-28]. The same is true for spin systems. When A C L, the sample is divided 
into correlated regions of size L c <C L, but such that kL c <C I. 

Inside of a correlated region, the average < S~ > weakly depends on the site j, so 
that this average can be taken out of the sum. Thus, we come to the definition of the 
collective radiation rate 

Tij) (26) 



(28) 

is the single-spin natural width. Inside such a correlated region, one has 

— — 2 < ft, ■ Hj > = (7 r — i5u)u . 

Averaging the spin equations of motion (10) over the spin degrees of freedom, we 
decouple the binary spin expressions < SfS? >, for % ^ j, as < S*f >< >, treating 
the local fluctuating fields (6) as random variables. At the same time, the spin products 
for the coinciding sites must be treated with caution. Averaging the last term in the first 
of Eq. (10), we employ the decoupling 

< srs; + s*s~ > = (2 - I) < s~ >< si > , 



= smifcr^j 



krij 



and of the collective frequency shift 



^ cos(A;r ii ) 

j M ij 

in which 

70-J- ^Sk* 
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which accurately interpolates the spin behaviour between the two exactly known limits 
of S = 1/2 and S -> oo [4,10,29]. 

Using inequality (25), the collective radiation rate (26) and the collective frequency 
shift (27) can be simplified. The retardation effects, due to the unit-step function entering 
Eqs. (26) and (27), manifest themselves at the time scale L c /c. Taking for L c the maximal 
value L c ~ 1/k, the retardation time can be approximated by Thus, Eqs. (26) and 
(27) can be represented as 

O AT 

lr - 7 oiV c (l - e-') , 5u * 7o ^ (l - e^) . (29) 

The time 1/ou is very short, usually it is the shortest among all other relaxation times. 
Therefore, for cut 3> 1, and taking into account that N c = pL\, we have for the collective 
radiation rate 

lr = loN c = J- p 2 Sk 3 N c (30) 



and for collective frequency shift 



8u -^t e = \ P»l S ( kL c? , (31) 

where p = N/V = N c /V c is spin density. Expressions (30) and (31) were, first, found 
by Ginzburg [30], who used a classical picture, which also was discussed in Refs. [4,31]. 
Here we have presented a fully quantum-mechanical derivation of these expressions and, 
in addition, have obtained their generalizations (26), (27), and (29). 

In studying spin dynamics, it is necessary to take account of the longitudinal and 
transverse relaxation rates. The former is due to spin-lattice interactions, and is denoted 
as 7i. The transverse relaxation rate is caused by the presence of dipole spin interactions 
and its standard value is 

72 = f PPIJS(STT) , (32) 

where n is the number of nearest neighbours. Quantity (32) is derived [19-21] under the 
assumption of a small longitudinal spin polarization. When the latter is large, so that 
the average polarization (24) substantially differs from zero, then the effective transverse 
relaxation rate narrows, becoming dependent on s. This effective narrowed rate was 
calculated and described in detail by Abragam and Goldman [21], who express its value 



72(s) 

through the moments 



CMl(s) 



M 4 (s) - Mf (s) 



1/2 



(33) 



M 2 (s) = M 2 (0)(1 - s 2 ) , M 4 (s) = 2.18M 2 2 (0)(1 - s 2 )(l - 0A2s 2 ) . 

The constant C in Eq. (33) depends on the line shape, being C = tt /2 for the Gaussian 
and C = 7r 2 /2 for Lorentzian lines. Substituting these moments in Eq. (33) and keeping 
in mind that s 2 < 1, we obtain 

72 ( S )=7 2 (0)(1- S 2 ), (34) 
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where 72(0) = yCM 2 (0) = 72. The total transverse relaxation rate, including the inho- 
mogeneous broadening 7^, is 

r 2 = 72(s)+7 2 *- (35) 

Averaging Eqs. (10) over the spin degrees of freedom, we take into account all de- 
scribed attenuations. Then, using the notation for the anisotropy frequency 

u D = \{2S-l)D (36) 

and for the effective spin frequency 

u s = u - lu d s , (37) 

we come to the evolution equations 

du 

— = -i{uj s + 4o - ^r 2 )w + fs , 

^ = -2F 2 w + (u*f + f*u)s , 
1 

■ d - t =-2 {u * f + ru) - ll{s - Ch (38) 
in which the effective force is 

/ = (7r - iSu)u + C • 

Because of the occurrence of the random variables £ an d £, corresponding to local spin 
fluctuations, these are stochastic differential equations. Averaging over the random fluc- 
tuations, treated as a Gaussian white noise [32], with the width 73, we follow the scale 
separation approach [5-10,14,22]. In this way, we obtain the equations for the guiding 
centers 

— = -2(72 + 72 - 12S 2 - lrS)w + 273S 2 , 

ds 

— = - lrW - 73 s - 7l ( s - C) , (39) 

containing all relaxation rates discussed above. 

Comparing the relaxation rates (30) and (32), we have 



72 3n V S + 1 

According to condition (25) for the existence of correlated regions, kL c 1, and since 
n ~ 10, the radiation rate 7 r is always much smaller than the dipole relaxation rate 72. 
As is clear from Eqs. (39), the value 7 r 72 plays no role in spin dynamics. Even taking 
the maximal correlation length L c ~ 1/k, one has 7 r /72 ~ 0.1, which does not noticeably 
influence the spin relaxation. Equations (39) contain no solutions corresponding to spin 
superradiance. 



10 



3 Numerical analysis and conclusion 



We analyzed numerically various solutions to Eqs. (39). The longitudinal relaxation 
rate was assumed to be small, 71 <C 72, which is usually the case. Other relaxation 
rates, expressed in units of 72, were varied in a wide diapason: 72 G [0,1], 73 G [0,1], 
and 7 r G [0,1]. Different initial conditions w = w(0) and s = s(0) were considered, 
satisfying the inequality Wo + s 2 , < 1. The solutions for 7 r = and / y r = 0.1 were 
found to be practically indistinguishable. Even unrealistically large 7 r = 1 resulted in 
insignificant changes, as compared to solutions with 7 r = 0. The main rates, governing 
the spin evolution, are 72, 72, and 73. In Figs. 1-4, we show the solutions for w(t) 
and s(t) at different relaxation rates, demonstrating the role of the latter, for the same 
initial conditions wo = and so — 1. These figures display the effect of self-induced 
dynamical coherence caused by the existence of dipole interactions. The main term in 
Eq. (39), producing the coherence intensity, is that containing 73. Without this term, 
no maximum in w appears. The term with 72 s 2 , due to the narrowing of the effective 
transverse relaxation, which results from a high spin polarization, is much less important 
and influences the spin motion by about 10%. If at the initial moment of time, an external 
transverse field imposes a noticeable coherence intensity Wq > 0, then we have the standard 
spin induction, as shown in Fig. 5. Both the spin induction and self-induced dynamical 
coherence occur on the time scale of the order of T 2 . Therefore, though these are coherent 
phenomena, they have nothing to do with superradiance which requires the pulse time to 
be much shorter than T 2 . Also, there is no spin reversal typical of superradiance. 

Recently, there have been attempts to detect the electromagnetic radiation generated 
in the avalanches of magnetization reversal in the molecular magnet Mn 12 -acetate, caused 
by the inversion of an external magnetic field [33]. However the most precise recent 
experiment [34] was unable to detect any significant radiation at well-defined frequencies. 
In any case, whether there is some electromagnetic radiation in the process of repolarizing 
Mni2 or it does not exist, it has nothing to do with superradiance. 

The accurate theory, based not on a cartoon model but on a microscopic consideration 
and taking into account the main interactions always existing in real molecular magnets, 
shows that the spin interactions through the common radiation field can never produce 
spin superradiance. This interaction is negligible as compared to dipole interactions and 
does not influence spin motion at all. Despite some analogy of spin systems with atomic 
ensembles displaying magnetodipole radiation [35] , there are principal differences between 
them. In atomic systems, all transverse relaxation rates are caused by the same electro- 
magnetic interaction, while in spin systems there ate several types of interactions, with 
the direct dipole interactions being the strongest and the most important. So, superradi- 
ance does exist in atomic systems. However it is unable to arise in spin systems, if these 
are not coupled to a resonator electric circuit. But if a molecular magnet is coupled to 
such a resonator, then spin superradiance can become achievable, as is suggested in Ref. 
[10] and described in detail in review [4]. 
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Figure Captions 



Fig. 1. Coherence intensity w(t) (solid line) and spin polarization s(t) (dashed line) 
as functions of time (measured in units of T 2 = I/72), with the initial conditions Wo = 
and s = 1, for the parameters 72 = and 73 = 0.1. 

Fig. 2. Coherence intensity (solid line) and spin polarization (dashed line) versus 
time (in units of T 2 ), with the initial conditions wo = and s = 1 5 for 7^ = and 73 = 1. 

Fig. 3. Coherence intensity (solid line) and spin polarization (dashed line) versus 
time (in units of T 2 ), with w = and s = 1, for 7 2 = 1 and 73 = 0.1. 

Fig. 4. Coherence intensity (solid line) and spin polarization (dashed line) versus 
time (in units of T 2 ), with Wo = and so = 1, for 7! = 1 and 73 = 1. 

Fig. 5. Spin induction. Coherence intensity (solid line) and spin polarization (dashed 
line) versus time (in units of T 2 ), for 72 = and 73 = 1, with the initial conditions 
wq = 0.75 and sq = 0.5. 
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